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FAST INTERPOLATION OF 
LARGE COLOR LOOKUP TABLES 

Background of the Invention 

The present invention relates to transforming a color from a first 
color space to a second color space. It finds particular application in conjunction 
with interpolating within a three-dimensional lookup table and will be described 
with particular reference thereto. It will be appreciated, however, that the 
invention is also amenable to other like applications. 

Color correction or transformation of color images is performed to 
convert images from one color space to another. Images may be represented using 
a device independent color space or a device dependent color space. Device 
independent color spaces represent color images independent from particular input 
and output devices so as to enable independent color imaging between input and 
output devices. Generally, a color image is transformed to a device dependent 
color space of an output device before the color image is rendered to ensure that 
the colorimetry of the original image is consistent with its reproduction. Such a 
transformation is performed regardless of whether the original image is defined in 
a device dependent color space or in a device independent color space. Closed 
systems that consist, for example, of a scanner and a printer that are not 
interconnected with external devices do not require an original color image to be 
transformed to a device independent color space. Such closed systems have color 
transformation systems that generate a color image represented by a first color 
space and subsequently convert that color space to a second color space before 
reproduction. Color transformation systems capable of independent color imaging 



represent color images in a third or a device independent color space which is 
readily transferable between a plurality of image output devices. 

Color images can be generated by an input image terminal such as a 
scanner or a color image creation program operating on a color workstation. Color 
5 images processed by a scanner or a workstation consists of a two-dimensional 
array of picture elements (pixels). The color of each pixel of an image may be 
represented using a plurality of color spaces. Scanner output is commonly 
transformed to a color space of tristimulus values, for example, additive primaries 
red, green and blue (RGB) color space. These values are typically a linear 

1 0 transformation of the standard XYZ coordinates of CIE color space, or a corrected 
transform of those values. In the-case of computer generated images, the computer 
workstation is operated by a user to create, edit, or view "softcopy" color images 
on the color monitor of the workstation. Colors selected by the user at the user 
interface of a workstation can also be defined in a color space of tristimulus values 

1 5 such as additive primaries RGB . 

An image generated by an image input terminal must be converted 
to subtractive primaries cyan, magenta, yellow and black (CMYK) or (simply the 
CMY) color space before being reproduced on an image output terminal such as a 
printer. CMY color space is typically used to represent the formulation of colored 

20 dyes, inks, or toners on paper. Printers typically operate by adding multiple layers 
of ink or colorant on each page. The response of the addition of colors by the 
printer tends to be relatively non-linear. Consequently, colors are defined for a 
particular printer and accordingly color spaces defined for a printer are device 
dependent. Thus, a printer receives information in a device independent color 

25 space from, for example, a workstation. The printer must then convert that 
information into its device dependent color space. There exist many different 
methods to convert between color spaces of images produced at a scanner or a 
workstation to a color space of images to be produced at a printer. 

Color correction consists of mapping points from a three- 

30 dimensional color space to a three-or-m ore-dimensional color space. This mapping 
depends on the physical properties of a marking device or printer system which is 
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typically nonlinear (as noted above). An effective approach to solving this 
problem is to use a coarse three-dimensional lookup table with interpolation. The 
lookup table provides an arbitrary mapping between different devices and the 
interpolation reduces the size of the table so that large amounts of memory are not 
5 required by the system to store a large number of sample points. In general, a 
lookup table contains values on a uniform grid for each of the three color 
coordinates of a color space. A common approach to interpolation is tri-linear 
interpolation which gives a linear weighting of all eight vertices. Alternatively, in 
tetrahedral interpolation the parallelepiped is divided into tetrahedra. The points 
10 within each tetrahedron are interpolated based on its four vertices. While 
tetrahedral interpolation is adequate for finding individual points at a lookup table 
of 16 3 nodes, it performs unnecessary computations when used with tables of 64 3 
nodes or more. 

More specifically, because color is defined with three (3) variables, 

15 a three-dimensional lookup table is commonly used. In RGB space, at a scanner or 
computer, space can be defined as three-dimensional with black at the origin of a 
three-dimensional coordinate system 0, 0, 0, and white at the maximum of a three 
dimensional coordinate system which in an 8-bit system, would be located at 255, 
255, 255. Each of the three (3) axes radiating from the origin point, therefore, 

20 define red, green, and blue, respectively. A similar construct can be made for the 
printer, with axes representing cyan, magenta, and yellow. In the 8-bit system 
suggested, there will be, however, over 16 million possible colors (i.e., 256 3 ), 
which are too many values for a 1:1 mapping of RGB to CMY within a reasonable 
amount of space. Accordingly, only a relatively small number of samples are 

25 conventionally used to do the mapping from RGB to CMY. 

The lookup table is conventionally constructed using slow methods 
that combine measured data with mathematical modeling. The time required to 
construct such a table of size 16 3 ranges from several minutes to one or two hours, 
depending on the method used. Constructing a large table of size 25 6 3 (where no 

30 interpolation is required) or even 64 3 , requires 4096 or 64 times as long, 
respectively. For example, if it takes 1 second to compute an entry, it will take 68 
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minutes to compute the 4096 entries of a 16 3 table. At that rate, it will take 194 
days to compute the 2 24 = 16,777,2 16 entries in a 256 1 table, and 72 hours, 49 
minutes to compute the 2 18 =262,144 entries in a 64 3 table. Thus, computing a table 
at such high density using conventional techniques is impractical. 
5 The present invention provides a new and improved method and 

apparatus which overcomes the above-referenced problems and others. 

Summary of the Invention 

A method for converting a specified color value from a first color 
space to a second color space identifies the specified color value in the first color 
10 space. A converted color space value is received from a final lookup table. The 
converted color space value is previously determined as a function of the specified 
color value and a mid-point interpolation and represents the specified color in the 
second color space. The converted color space value is stored in a memory device. 

In accordance with one aspect of the invention, before the converted 
15 color space value is received, the final lookup table is created as a function of an 
initial lookup table, which has a lower resolution relative to the final lookup table 
and includes a plurality of the conversion values. 

In accordance with another aspect of the invention, a number of 
mid-point interpolations is performed for a current color of interest as a function of 
20 the resolution of the final lookup table. 

In accordance with a more limited aspect of the invention, the mid- 
point interpolation is performed by determining an anchor vertex, which is defined 
within the final lookup table, as a function of the current color, and determining a 
current vertex, which is defined within the final lookup table, as a function of the 
25 current color and the anchor vertex. A current intermediate color is determined by 
performing a mid-point interpolation as a function of the anchor and current 
vertices. 

In accordance with an even more limited aspect of the invention, if 
more than one mid-point interpolation is performed, a second current vertex, which 
30 is defined within the final lookup table, is determined as a function of the current 
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color, the anchor vertex, and the current vertex. A second current intermediate 
color is determined by performing a mid-point interpolation as a function of the 
second current vertex and one of the anchor and first current vertices. 

In accordance with another aspect of the invention, the anchor 
5 vertex and the current vertex are determined by determining a vertex, which is 
defined within the final lookup table, as a function of a node in the first color space 
having largest coordinates less than corresponding coordinates of the specified 
color. 

In accordance with another aspect of the invention, if the final 
10 lookup table is less than a predefined resolution, a plurality of intermediate color 
space conversion values are received from the final lookup table. The converted 
color space value is determined as a function of the intermediate color space 
conversion values. 

In accordance with a more limited aspect of the invention, if the 
15 predefined resolution is half the resolution of the input colors, two conversion 
values are received from the final lookup table. 

An advantage of the present invention is that it provides a relatively 
fast technique for finding individual points at a granularity of 256x256x256 within 
a lookup table having 64x64x64 nodes. 
20 Still further advantages of the present invention will become 

apparent to those of ordinary skill in the art upon reading and understanding the 
following detailed description of the preferred embodiments. 

Brief Description of the Drawings 

The invention may take form in various components and 
25 arrangements of components, and in various steps and arrangements of steps. The 
drawings are only for purposes of illustrating a preferred embodiment and are not 
to be construed as limiting the invention. 

FIGURE 1 illustrates a system for converting a color value from a 
first color space to a second color space according to the present invention; 
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FIGURE 2 illustrates a flowchart for converting a color value from 
a first color space to a second color space according to the present invention; 

FIGURE 3 illustrates a flowchart for determining a lookup table 
according to the present invention; and 
5 FIGURE 4 illustrates a portion of a color space according to the 

present invention. 



Detailed Description of the Preferred Embodiments 

To simplify discussion in the following, we shall refer to number of 
inter-node spaces along one dimension of a lookup table as the table resolution. 

10 Thus a 17x17x17 table, which has 17 3 = 4913 nodes, has a resolution of 16. The 
table may also be viewed as containing 16x16x16 cubes, each containing 8 nodes. 
All of the cubes except those at the gamut boundaries share all of their vertices 
(nodes) with neighboring cubes. Those at the gamut boundaries share some but not 
all of their vertices with neighboring cubes. A given cube is identified by the 

15 location of its lower left, front vertex (i.e. the vertex with the lowest coordinate 
values). Indices into the table for all nodes that are the lower left, front vertex for 
some cube require 4 bits, so, the table is said to have 4 bit resolution. 

With reference to FIGURES 1 and 2, a system 10 for converting a 
specified color value (pixel value) from a first (initial) color space to a second 

20 (final) color space includes an input device 12 (e.g., a scanning device), a display 
device 14, a memory device 16, a processing device 20, and an output device 22. 
The processing device 20 communicates with both the memory device 16 and the 
output device 22. In the preferred embodiment, the memory device 16 is included 
within the processing device 20; however, other configurations are also 

25 contemplated. It is to be understood that the first and second color spaces may be, 
for example, device dependent color spaces (e.g., RGB or CMYK) or colorimetric 
color spaces (e.g., L*a*b*). Furthermore, the display device 14 may, for example, 
operate in the RGB color space, while the input and output devices 12, 22, 
respectively, operate in the CMYK color space. 
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A method 24 for converting the specified color value (pixel value) 
from the first color space to the second color space begins in a step A. A three- 
dimensional lookup table, preferably stored in the memory device 16, for 
converting the color value between the color spaces is assumed to have three or 
5 four bit resolution (i.e., 9x9x9 or 17x17x17 nodes). As will be described in more 
detail below, a final lookup table, which has six, seven or eight bit resolution (i.e., 
65x65x65, 129x129x129 or 257xx257xx257) nodes, is determined in a step B. 
The final lookup table is stored in the memory device 16 in a step C. Because the 
final lookup table has more nodes than the initial lookup table, the final lookup 

10 table has a higher resolution relative to the initial lookup table. The nodes of the 
lookup tables include conversion values (transformation values) for converting 
(transforming) respective color values (pixel values) from the initial color space to 
the final color space. 

Once the final lookup table is determined, an image 30 is scanned, 

15 in a step D, into the memory device 16. It is to be understood that the image 30 is 
scanned in a first color space (e.g., a device dependent color space such as RGB) 
and converted to a second color space (e.g., a colorimetric color space such as 
L*a*b*) before it is stored in the memory device 16. Furthermore, the image data 
is converted on a pixel-by-pixel basis. 

20 A current pixel value (specified color) to be converted is identified, 

in a step E, within the image 30. A conversion value (transformation value), which 
is selected as a function of the current pixel value, is received into the processing 
device 20 from the final lookup table in a step F. If the resolution of the final 
lookup table is less than eight bit, multiple intermediate conversion values are 

25 selected, and processing continues with a step G; otherwise a single conversion 
value is selected, and processing passes to a step H. In step G, a single conversion 
value is calculated from the multiple conversion values. The single conversion 
value is stored in the memory device 16 in step H. In this way, the current pixel 
value is converted from the first color space to the second color space, and stored 

30 in the memory device 16. 
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It is to be understood that a similar conversion process may be 
implemented for converting the current pixel value to the appropriate color space 
for the output device 22. In this situation, the current pixel value is optionally 
produced, in a step I, via the output device 22. The output device 22 is preferably 
a digital output device (e.g., printing device), which operates in a xerographic 
environment, which produces either black-&-white or color output. 

A determination is made, in a step J, whether more pixel values 
remain to be converted in the image 30. If more pixel values are to be converted, 
control returns to the step E for identifying the next current pixel value. If more 
pixel values are not to be converted, control passes to a step K for stopping the 
process. 

With reference to FIGURE 3, the step G of converting the current 
pixel starts in a step Gl. In one preferred embodiment, the resolution of the final 
lookup table is assumed to be 2 7 (129x129x129 nodes). A current color of interest 
is identified in a step G2. 

With reference to FIGURES 3 and 4, the final lookup table is 
represented as a cube 40 of 129x129x129 nodes (vertices). Then, the only 
locations (vertices) of concern (other than at the gamut boundaries) are those points 
42a, 42b, 42c, 42d, 42e, 42f, 42g, 42h (vertices) marked with spheres. Other 
points are shared with other like cubes. It is clear from FIGURE 4 that, if the final 
lookup table has 7 bit resolution, then to find any point at 8 bit resolution, it is 
sufficient to average one of the vertices 42a, 42b, 42c, 42d, 42e, 42f, 42g, 42h (an 
anchor vertex) with one of the anchor vertex's seven (7) neighbors. For example, 
the lower, left vertex 42a is chosen as the anchor vertex in a step G3. Then, one of 
the seven (7) neighboring vertices 42b, 42c, 42d, 42e, 42f, 42g, 42h, is chosen, in a 
step G4, as a current vertex (e.g., 42b). The anchor vertex 42a is the closest node 
with coordinates in the initial color space less than those of the current color of 
interest. The current vertex 42b is selected so that the current color of interest lies 
midway along a line from the anchor vertex and the current vertex 42b. That is, it 
is the node that identifies the cube containing the current color of interest. An 
intermediate vertex (intermediate color) 42i is determined, in a step G5, via a mid- 



point interpolation between the anchor vertex 42a and the current vertex 42b. The 
process stops in a step G6. 

In a second embodiment the final table resolution is 2 6 , or 6 bits. A 
simple method of determining the conversion is to first locate a pair of 
5 intermediate nodes. These intermediate nodes are those nodes in a hypothetical 2 7 
bit grid that would contribute to the computation of the converted pixel value for 
the current color of interest. These nodes are then calculated in a method 
analogous to the method for converting from 7 to 8 bit resolution. To obtain the 
final color, the intermediate nodes are used to convert from 7 to 8 bit resolution. 

10 This simple conversion apparently involves three times the operations that are 
required to go from 7 to 8 bit resolution: two times to get the two intermediate 
nodes, and once to use them. However, as will be shown below, it is possible to 
reduce the cost to less than twice the cost of going from 7 to 8 bit resolution. 

The mid-point interpolation averaging operations discussed in the 

15 present invention involve significantly smaller computational cost than tetrahedral 
and other conventional interpolation methods. A mid-point p between two points q 
and r is given mathematically by p = (q + r) / 2. In a software implementation, this 
may be represented as p= ((q+ r+l)» 1). Hence, a midpoint subdivision 
requires two (2) adds and a shift. If we wish to compute two (2) mid-point 

20 averaging operations, we have p= (s+ (q+ r)/ 2)1 2), and this translates in 
software to p = (q + r + (s « 1) +2) » 2) (which requires three (3) adds and two 
(2) shifts). Furthermore, such a computation is more accurate on average than two 
(2) separate mid-point subdivisions. Note that mid-point averaging requires no 
multiplications, a significant gain over conventional interpolation techniques. As 

25 we will see below, the 3-D interpolation for each input color may be performed 
with no more than two (2) mid-point averaging operations per output separation. 

As discussed above, the cube shown in FIGURE 4 shows a 
64x64x64 tiling of color space (that is, there are 64 cubes, or 65 nodes in each 
direction). Points anywhere in a 7 bit representation are determined via mid-point 

30 interpolations between two (2) nodes within the 64x64x64 tiling. Of all points 
having eight bit representations, those with the low order bit 0 lie both on the 
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257x257x257 grid, and on the 129x129x129 grid. These represent approximately 
one in every 8 points. Points with 8 bit resolution anywhere else are determined 
via mid-point interpolations between two (2) nodes within a hypothetical 7 bit 
resolution representation, at least one of which is a vertex in the 64x64x64 tiling. 
5 Because at least one vertex of the original 64x64x64 tiling can be used in the final 
mid-point interpolation and the other vertex used in the final mid-point 
interpolation is calculated using mid-point interpolation from the original 
64x64x64 tiling, the method requires only two mid-point interpolations. In this 
manner, there are always at most two (2) averagings, which involve at most three 

10 (3) points from the 64x64x64 tiling. 

In one case out of 64, no averagings are needed, but to simplify 
logic, this case is done by averaging the point with itself. In that and an additional 
seven (7) cases, the result of the first averaging yields the desired result. 

The implementation of the mid-point averaging technique assumes a 

15 final color lookup table having a six bit resolution is given or has been computed. 
A single point in, for example, the RGB color space is interpolated by using the six 
(6) high order bits to obtain cell indices j, k, and 1. The low order bits are 
concatenated according to Lowlndex = (R & 0x3) | ((G & 0x3) « 2) | (B « 4). 

The selection of which nodes are used in the averaging is not easily 

20 suited to rapid computation. For example, all the nodes represented as spheres in 
FIGURE 4 may be computed by averaging two nodes at vertices of the cube. 
While it is not difficult to compute which node should be the other node in the 
average, once node 42a is selected as an anchor node, it is faster to store that 
information in a lookup table, indexed on the low order bits of the node to be 

25 computed. As a second example, consider the point (not shown), midway between 
nodes 42c (on the front left edge) and 42e (along the main diagonal). That node 
will be computed as the average of node 42a and an intermediate node (not shown) 
at the center of the top face of the cube. It is even more attractive to use lookup 
techniques in order to select these three nodes and their order of interpolation than 

30 for the first example. 
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In light of the foregoing, all the bits of Lowlndex are used as an 
index into a 256 element lookup table. Entries in the 256 element lookup table 
repeat every 64 nodes. Ten (10) one (1) bit numbers bO ... b9 (which may be stored 
as bytes in a software environment) are retrieved from the table. If bO is 0, then 
5 only one average is needed. The required average would then be retrieved from the 
lookup table as a function of ((LUT[j + bl] [k + b2] [1 + b3] + LUT[j + b4] [k + 
b5J [1 + b6] + 1) » 1 (7/64ths of the time), where LUT represents "lookup table." 
Otherwise, two (2) averaging operations are required (q + r + (s « 1) +2) » 2): 
((LUT[j + bl] [k + b2] [1 + b3] + LUTJj + b4] [k + b5] [1 + b6] + LUT[j + b7] [k + 

10 b8] [1 + b9] «1) + 2 » 2. This version involves one (1) branch (depending on 
whether two (2) or one (1) interpolations are needed), and so is amenable to 
machines for which branching is expensive (as is typical of modern machines). An 
alternative, with still one (1) branch, but this time through a jump table, and one (1) 
fewer lookup, is to use Lowlndex as the control variable for a switch statement. In 

1 5 summary, the operation counts required are given below, where N is the number of 
output dimensions: 



Operation 


Tetrahedral 
Interpolation 


Method 1 


Method 2 


Add 


3N 


(2 57/64)N 


(2 57/64)N 


Shift 


N 


(1 57/64)N 


(1 57/64)N 


Multiply 


3N 


0 


0 


3D Lookup 


4 


3 


3 


ID Lookup 


0 


1 


0 


N-way Branch 








(Switch) 


0 


0 


1 


Branch 


2.5 




0 



It is anticipated that the savings in branches, multiplies, and three- 
20 dimensional lookups will easily outweigh the extra shifts and the switch or one- 
dimensional lookup. 

Note that while the technique described above requires a 64x64x64 
lookup table, it is not necessary that such large lookup tables be persistent in a 
color management system. One can easily upsample lookup tables of conventional 
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sizes (e.g., 16x16x16 to 64x64x64) just prior to processing the image. In fact, the 
technique just described above provides one (1) possible way to repeatedly 
upsample the lookup table by a factor of two (2) in each dimension until the 
desired sized is reached. However, the method of interpolation, because it always 
5 uses a limited number of vertices from a cube, is unlikely to provide satisfactory 
results when applied recursively. Therefore, a method is provided for upsampling 
a table using linear interpolation or cubic interpolation. 

It is to be understood that a sparse table is constructed using 
traditional (data- or model-based) techniques; the full table is then built using 

10 interpolation. The interpolation uses a fast bisection technique. Either tri-linear or 
tri-cubic interpolation may be used depending on the application. Typical table 
sizes are 9x9x9 or 17x17x17. While the method as described above uses fast 
interpolation on the final table, the interpolation using tetrahedral or trilinear 
interpolation could also be used, at a sacrifice of speed. 

15 As discussed above, constructing high density 3-4 color conversion 

tables is slow (a typical characterization, yielding a 17x17x17 table using current 
techniques on a typical workstation may take nine (9) minutes, but with multiple 
rendering intents, print modes, media, etc. the full characterization may run into a 
few hours). The characterization is relatively smooth, meaning with a cubic 

20 interpolant an 8x8x8 table is likely sufficient for even the most demanding 
applications. It is even possible for a high end application to generate a full table 
(i.e. 256x256x256 — 2 24 entries), but that would take time proportional to the 
number of entries, and if a 17 3 table takes 9 minutes, (roughly 2 12 entries), then a 
full table would require 4000 times as long, or the better part of a month. 

25 Small lookup tables (e.g. 5x5x5) exhibit two problems: first, they 

may fail to capture the inherent non-linearity of the conversion process (this will 
depend on the device, manifesting itself as large errors in regions where the 
underlying function has high curvature); second, they may show artifacts due to the 
linear interpolation occurring between nodes (potentially generating artifacts from 

30 kinks in the interpolant making sweeps less smooth). Moderate sized tables (e.g. 
9x9x9) may continue to have linear interpolation artifacts, but better capture device 



behavior. A table at 1 7x 1 7x 1 7 contains all of the information content of a typical 
characterization, however it is possible that linear interpolation may still be a 
problem for the most demanding applications. Beyond 33x33x33, linear 
interpolation is highly unlikely to cause image quality problems. However, very 
5 large tables reduce or eliminate the need for interpolation. A 256x256x256 table 
converting from RGB to L*a*b* (for a scanner), returns the desired color in a 
single table lookup, but requires 48 Mbytes of memory. A 129x129x129 table 
requires only 6 Mbytes of memory (and two adds and a shift to interpolate). 
Because each interpolation is a midpoint subdivision, substantial computational 

1 0 savings is had even with the 1 29x 1 29x 1 29 table. 

Because of the value of using large tables, interpolated or not, the 
present invention addresses constructing such large tables. Midpoint subdivision is 
a known, fast means of computing curves and surfaces; this is a generalization to 
volumes. Midpoint subdivision is numerically stable, but not the only way of 

1 5 doing the interpolation. For example, one could use forward differencing to build a 
table with any interpolant to any size; for large scale factors this would be faster. 
However, midpoint subdivision is fast enough for most applications, as timing tests 
(described below) have shown. 

20 Linear Interpolation 

If the input table is at high enough density that the only 
improvement from going to the large table is speed, the conversion from the input 
table size to the final table size may be done by repeated midpoint subdivision as 
follows: 

25 Assume the table to be computed (the final table) has 2" + 1 entries 

per dimension, and the input table has 2 m + 1 entries, n>m. In a floating point 
implementation, begin by storing the table sparsely, so that each entry of the input 
table is in its final location. Then: 
For each dimension 

30 Fork = mton-l 
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For each pair of existing entries e^ e i + 2 n ' k ; 

Form a new entry e ( + 2 n - M = 0.5 * (e-, + e, + 2 "- k ) 
Here i is the index in the current dimension, k is the current 
logarithm (base 2) of the table size; and other dimensions are held fixed. 

At each iteration of the For k loop, the table size doubles in the 
current dimension; this continues until the table reaches its final size in the current 
dimension. This is repeated for each dimension. Thus the table might grow from 
9x9x9 to 9x9x17 to 9x9x33 to ... to 9x9x257, then to 9x17x257, ...and finally 
from 129x257x257 to 257x257x257. 

In fixed point, the values should be pre-loaded with n-m guard bits. 
Then the expression for forming a new entry is simply: 

e i + 2 n - k -' = (e i + e i + 2 n - k )»l. 
Once the table is fully computed, the guard bits must be discarded: 
e s + 2 n - ra -' = (e ; + 2 n - m - , +2 n - m - 1 )»(n-m); 

A test implementation found that going from 17x17x17 to 
257x257x257 required only 13 seconds per separation. 

Cubic Interpolation 

For smaller input tables or higher accuracy, cubic interpolation may 
be used. The Bezier cubic provides interpolation of the nodes, and with properly 
arranged "intermediate nodes", CI continuity. Other methods could also be used. 
In one dimension, the process proceeds in two steps: 

1 . Generate intermediate nodes as: 
n 4i = Pi 

n 4i+ i =Pi + (p i + i.Pi-i/4 
n 4l+ 2 = Pi + . +(Pi-P i+2 / 4 
n 4i+ 3 = P, + i 

2. For each level 

For i = 0, 3, 6, ... subdivide the segment nj..n i+3 , as follows: 
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Compute nj i+ , as (nj+n i+ |)/2, and similarly compute 
nj+1,,+2 and iv 2 i+3 . 

Average adjacent pairs from the three just computed 
- yielding i+2 and rv u+3 ; 
5 Average the pair just computed, yielding iij i+3 . 

Replace the segment with the two segments ri;, n^,, 
^,1+2, ni, i+3 ; and n^, n; + , ji+3 , rv 2 i+3 , n i+3 



At the last iteration, only the new segment joint need be retained. 
10 This would be repeated in each dimension as for trilinear 

interpolation. 3 guard bits per level of subdivision are needed for accuracy here, in 
a fixed point implementation. 

This is approximately six times as expensive as linear interpolation, 
which is still very fast (roughly 80 seconds per separation for a full 257x257x257 
15 table.) 



Applications 

One application already mentioned is the very large table 
(257x257x257 or 129x129x129). Another application is fast recharacterization of 

20 a device. If a device is well characterized by the information in an 8x8x8 table, but 
the system uses 16x16x16 or even 32x32x32 to avoid linear interpolation artifacts, 
the 8x8x8 table can be computed in 1/8 the time of a 16x16x16, or 1/64 the time of 
a 32x32x32 (in other words a few minutes), and then scaled up with cubic 
midpoint interpolation in a fraction of a second (even for the 32x32x32 case). 

25 The invention has been described with reference to the preferred 

embodiment. Obviously, modifications and alterations will occur to others upon 
reading and understanding the preceding detailed description. It is intended that 
the invention be construed as including all such modifications and alterations 
insofar as they come within the scope of the appended claims or the equivalents 

30 thereof. 
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